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The deformation of a nonrotating body resulting from the application of a tidal field is measured 
by two sets of Love numbers associated with the gravitoelectric and gravitomagnetic pieces of the 
tidal field, respectively. The gravitomagnetic Love numbers were previously computed for fluid 
bodies, under the assumption that the fluid is in a strict hydrostatic equilibrium that requires the 
complete absence of internal motions. A more realistic configuration, however, is an irrotational state 
that establishes, in the course of time, internal motions driven by the gravitomagnetic interaction. 

We recompute the gravitomagnetic Love numbers for this irrotational state, and show that they 
are dramatically different from those associated with the strict hydrostatic equilibrium: While the 
Love numbers are positive in the case of strict hydrostatic equilibrium, they are negative in the 
irrotational state. Our computations are carried out in the context of perturbation theory in full 
general relativity, and in a post-Newtonian approximation that reproduces the behavior of the Love 
numbers when the body’s compactness is small. 


I. INTRODUCTION AND SUMMARY 

A body subjected to an applied tidal field suffers a deformation that depends on the details of its internal structure. 
When the body is nonrotating, these details are encapsulated in a set of gravitational Love numbers Kf and 
and a measurement of the tidal properties of a body can reveal, through the Love numbers, important information 
regarding this internal structure. This observation has motivated the development of a relativistic theory of tidal 
deformation and dynamics, in the context of the measurement of tidal effects in gravitational waves emitted by 
neutron-star binaries [THU] and during the capture of solar-mass compact bodies by supermassive black holes [T5H^ . 
Tidal invariants have been incorporated in point-particle actions to account for the tidal response of an extended 
body [ 2 TH 25 ] . 

While the response of a self-gravitating body to an applied gravitoelectric tidal field is familiar from Newtonian 
theory (see, for example, Sec. 2.5 of Ref. |26| for a thorough treatment), its response to a gravitomagnetic tidal field 
is a relativistic effect that has no analogue in Newtonian gravity. This effect was first explored by Favata m in the 
context of post-Newtonian theory, and subsequently by Damour and Nagar [55] and Binnington and Poisson [55] in 
full general relativity. We examine it further in this work, and inspired by Favata, we lift an important restriction on 
the types of fluid configurations that were allowed in the earlier, fully relativistic work. 

The gravitomagnetic Love numbers of a fluid body were computed by Damour and Nagar and Binnington 

and Poisson under the assumption that the tidal interaction is sufficiently slow that it never takes the body out of 
hydrostatic equilibrium. This is a good approximation for many circumstances; for example, it is expected to hold for 
most of the orbital evolution of a compact binary system, up to the point where merger is about to take place. But the 
hydrostatic equilibrium considered in the earlier work is a strict one that forbids the existence of fluid motions within 
the body; the compact body is assumed to be strictly static, except for the parametric time dependence communicated 
by the slowly changing state of the tidal environment. 

Our main purpose in this paper is to point out that the strict hydrostatic equilibrium is too severe a restriction 
on the body’s internal physics. We follow instead Shapiro [3D] and Favata m, and take the fluid to be in an 
irrotational state that permits internal motions driven by the gravitomagnetic interaction with the tidal environment. 
We recalculate the Love numbers AT™*^® for this configuration, and show that they are dramatically different from 
those associated with the strict hydrostatic equilibrium: While the Love numbers are positive in the case of strict 
hydrostatic equilibrium, they are negative in the irrotational state. 

In our work the tidal field is still taken to vary slowly, and the fluid is still taken to be in an approximate hydrostatic 
equilibrium, in the sense that the fluid’s physical variables, such as density, pressure, and velocity field, carry only 
a parametric dependence upon time that reflects the slow evolution of the tidal environment. But internal motions 
are now allowed. As Shapiro and Favata have shown, these internal motions are a consequence of the conservation 
of relativistic circulation within the fluid, and are established whenever the tidal field exhibits a time dependence, 
however slow it may happen to be. On the other hand, the strict hydrostatic equilibrium adopted in the earlier 
work requires the tidal environment to be strictly stationary; it is a far less realistic description of the fluid. Our 
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considerations in this paper are limited to the gravitomagnetic interaction; as we shall show, the switch from strict 
hydrostatic equilibrium to the irrotational state has no impact on the body’s gravitoelectric response. 

We begin our developments in Sec. [IT] with a description of the unperturbed state of an isolated, self-gravitating 
body consisting of a perfect fluid; we take the unperturbed configuration to be static and spherically symmetric. In 
Sec. |III| we introduce a perturbation and examine the relativistic Euler equation that governs the perturbed state 
of the fluid. We continue the discussion in Sec. |IV] by working out the consequences of the relativistic circulation 
theorem for our perturbed configuration; we show that the irrotational state comes with internal motions that are 
forbidden in the strict hydrostatic equilibrium. 

In Sec. |V] we specialize the perturbation to a gravitomagnetic tidal field, and we calculate the body’s response to 
this field when the fluid configuration is in the irrotational state. The tidal environment is generic and characterized 
by an £-pole gravitomagnetic moment that is assumed to vary slowly with time. This Cartesian tensor is 

symmetric and tracefree (STF), and in a quasi-Lorentzian frame (t,x^) it appears in the time-space components of 
the metric tensor. 
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Here, Cjpq is the completely antisymmetric permutation symbol, is the gravitomagnetic Love number of degree 

£ > 2, M is the body’s gravitational mass, := SjkX^x’^, and dots indicate relativistic corrections of order 2M/r and 
higher; we work in relativistic units with c = G = 1. This expression for applies to a domain R < r < rout ^ b 
bounded internally by the body’s radius R and externally by an outer radius rout required to be much smaller than b, 
the distance to the external matter responsible for the tidal field. The first term in gtj, which grows as r^, represents 
the external tidal field, and the second term, which decays as represents the body’s response to the applied 

tidal field, quantified by The tidal moments Bkik 2 -ki{t) can be thought of as a collection of functions of time 

that cannot be determined by the Einstein field equations restricted to the domain R < r < rout; or they can be 
viewed as components of the Riemann tensor differentiated £ — 2 times and evaluated in the regime r ^ M, 
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in which the STF label instructs us to symmetrize the kik 2 ■ • • kg, indices and remove all traces. This notation, and 
our expression for the external piece of goj, is imported from Zhang’s pioneering work [31) . For a generic tidal 
environment the dominant moment is Bjk-, and the most relevant Love number is but a formulation of the 

body’s tidal response can be provided for any multipole order. 

The gravitomagnetic Love numbers depend on the details of the body’s internal structure, as determined by 

its equation of state, which we take to be of the zero-temperature form p = p{p), e = e(p), where p is the rest-mass 
density, p the pressure, and e the density of internal energy. For concreteness we adopt a simple pol^ropic model 
p = RTpi+i/", e = np, where K and n are constants. A sample of our computations is displayed in Fig. Id which plots 


as a function of 2M/R for £ = 2 and selected values of the polytropic index n. We observe that the tidal response 
of an irrotational fluid is dramatically different from the response of a fluid in strict hydrostatic equilibrium: the 
gravitomagnetic Love numbers of an irrotational body are negative, while they are positive for a body in strict 
hydrostatic equilibrium. Both sets of Love numbers, however, share the properties that they decrease (in absolute 
value) with increasing 2M/R and with increasing n; because increasing n decreases p for a given p, and therefore 
produces a body that is more centrally dense, both properties are associated with the fact that a more compact body 
develops sma ller multipole moments. 

we exploit the methods of post-Newtonian theory to calculate in the limit 2M/R —>■ 0, for any 
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equation of state. We obtain the simple expression 
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which can be evaluated explicitly once the equation of state is specified. The parameter A tracks the internal motions 
associated with the fluid’s irrotational state; setting A = 1 places the body in the irrotational state and gives rise 
to a negative for any £ and any equation of state, while setting A = 0 places the body in a strict hydrostatic 
equilibrium and gives rise to a positive gravitomagnetic Love number. The post-Newtonian values appear as circled 
data points in Fig.[^ and we see that they accurately reproduce the 2M/R -+ 0 limit of the relativistic curves. 
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FIG. 1. Rescaled gravitomagnetic Love number for 1 = 2. The Love numbers are computed for selected polytropes with 
polytropic index n, and plotted as a function of 2M/R, up to (and slightly beyond when n > 1) the maximum value for the 
given equation of state. The positive values correspond to the strict hydrostatic equilibrium. The negative values correspond 
to the irrotational state. The circled data points at 2M/R = 0 are the post-Newtonian values displayed in Eq. (1.41. 


We conclude this introduction with a restatement of our main message. A nonrotating body in a tidal interaction 
with remote matter is expected to be in an irrotational state for which internal motions get established over time. 
These internal motions have a dramatic influence on the body’s gravitomagnetic Love numbers: while they are positive 
in the strict hydrostatic equilibria examined in previous works, they are negative in the irrotational state. This effect 
is captured by the post-Newtonian expression of Eq. (1.4), which is valid in the limit 2M/R 0 for any equation of 
state. 


II. UNPERTURBED CONFIGURATION 

The unperturbed body is taken to be static and spherically symmetric, and to consist of a perfect fluid with rest- 
mass density p, pressure p, density of internal energy e, and density of total energy p = p + e. The background 
spacetime has a metric given by 

ds'^ = dt^ + f~^ dr'^ + dil'^ ( 2 . 1 ) 

with / := 1 — 2m/r, in which ^ and m depend on the radial coordinate r, and d^^^ := dO"^ + sin^ 9 d(j)^. The metric is 
a solution to the Einstein field equations with energy-momentum tensor 

T°‘P = {p + p)u°‘uf^ + pg°^^, (2.2) 

in which is the fluid’s velocity field, with only nonvanishing component u* = e~'^. The metric functions are 
determined by 

m + Attv^p 

j.2j: 


m! = 47rr^/i, 


(2.3) 
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in which a prime indicates differentiation with respect to r. In the vacuum exterior the field equations produce 
the Schwarzschild solution m = M = constant and = 1 — 2M/r. The body’s surface is situated at r = i?, as 

determined by the condition p{R) = 0. 

The conservation equation = 0 gives rise to the relativistic Euler equation, 

(/r + p)aa + {Sj + UaU^)V pp = 0, (2.4) 

where a“ := pu°‘ is the covariant acceleration, and the first law of thermodynamics. 


+ (m + p)^au‘^ = 0 . 


(2.5) 


For the unperturbed configuration the only nonvanishing component of the acceleration vector is 
reduces to 


and Eq. (2.4) 


p = -{p+p)ip' 


(/i + p){m + Anr^p) 


( 2 . 6 ) 


the Tolman-Oppenheimer-Volkov equation. 


III. PERTURBED CONFIGURATION 


We next allow the fluid configuration to be perturbed by a time-dependent, external tidal field. The metric becomes 
Pap+Pap, and the fluid quantities are shifted to p + Sp, p + Sp, u'^ + 6u°‘, Ua + Sua, and aQ-|-i5aQ,. We work consistently 
to first order in the perturbation, and note that Sua = Pap^u^ + Papu^■ 

Normalization of the perturbed velocity vector in the perturbed metric implies that UaSu°‘ = —^papu°‘u^^ or 
6u* = \e~^'^ptt- The remaining components are denoted Su^ := u’' and Su'^ := with 9"^ = {9,(j)). We let 
Va '■= Papv^ and find that the components of 6ua are given by 


6ut = ^e ^Ptt, 5Ur = Vr + e "^Ptr, 5uA = VA + e '^PtA- 
A straightforward computation further reveals that 

Sat = 

Sar = e~'>’dt{vr + e~*ptr) - ^dr{e~'^'^ptt), 

5aA = e~'f'dt{vA + e~'>'ptA) - 

The perturbed configuration is governed by Euler’s equation, which becomes 

{p -I- p)daa + {Sp + Sp)aa -f {S£ -I- UaU ^)VpSp + (UaSu^ + U^SUa)VpP = 0 


(3.1) 


(3.2a) 

(3.2b) 

(3.2c) 


(3.3) 


after the perturbation. With the information provided above, we find that the radial component reads 
e~'^{p+p)dt{vr + e~'*ptr) - ^{p + p)dr{e~'^'^ptt) + {Sp + 6p)'tjj' + drSp = 0, 


(3.4) 


while the angular components take the form of 


e '^{p+ p)dt {vA + e "^ptA) + Qa 


1 

2 


{p + p)e ‘^^ptt + Sp 


= 0 . 


(3.5) 


The time component of Euler’s equation returns a trivial 0 = 0. 


IV. IRROTATIONAL CONFIGURATION 

The circulation of a relativistic fluid around a closed curve c is defined by 


C(c) := 


hua da:“, 


(4.1) 
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where h is the specific enthalpy defined by d\nh := (/z + dp, and dx'^ is the coordinate increment along c. It is 
known that C is the same for any circuit c that surrounds a given fluid world tube. (The world tube can be thought 
of as a bundle of streamlines, defined as the world lines of fluid elements.) Thus, if Ci surrounds the world tube at 
a time t = ti, and if C 2 surrounds the same world tube at a time t = t 2 , then C(c 2 ) = C'(ci) and the circulation is 
conserved. A proof of the circulation theorem can be found in Synge’s 1937 review of relativistic hydrodynamics |32) . 

Because Ua dx°‘ = 0 for the unperturbed configuration, we have that the unperturbed fluid is irrotational: C(c) = 0 
for any purely spatial circuit c. The circulation of the perturbed configuration is then 

C(c) = j)[6hua + hSua)dx°‘ = <j> hSua dx°'. (4.2) 

We assume that the fluid begins in an unperturbed state, so that the perturbation vanishes at some initial time 
t = to- This implies that C(co) =0 for any circuit cq tangent to the hypersurface t = to. Conservation of circulation 
then guarantees that the perturbed fluid is irrotational at all times: C'(ct) = 0 for any circuit Ct tangent to any 
hypersurface t = constant > tg- And because the integral of h6uadx°‘ must vanish for any circuit Cj, we conclude 
that an irrotational fluid configuration must satisfy 


Sur = 0 = 6ua 


at all times. Importing Eq. (3.11, we have that 

Vr + e~'*Ptr =0, VA+ e~'^ptA = 0 


(4.3) 


(4.4) 


for an irrotational configuration. 


We next insert the second of Eqs. (4.4) within Eq. (3.5). The equation integrates to 




and this can be substituted back into Eq. (3.4), along with the first of Eqs. (4.4). This yields 




With the help of Eq. (2.6) these expressions become 


Sp = —rp'F, 6p = —rp'F, 


F := 


Ptt 


2re^’^'0' 


(4.5) 


(4.6) 


(4.7) 


These equations imply that a spherical surface of constant p or p at radius r in the unperturbed configuration 
is deformed by the perturbation to a nonspherical surface at r(l + F). Equations (4.7) are a consequence of the 
perturbed Euler equation and the conservation of circulation for an irrotational fluid. They allow the perturbation 
Pap to be time-dependent, but they are identical in form to the equations of hydrostatic equilibrium derived, for 
example, by Landry and Poisson [55] . 

We are interested in the deformation of a fluid body placed in a tidal environment, and we shall henceforth assume 
that the tidal field varies slowly with time. In this context, the state of the fluid is at all times an approximate 


hydrostatic equilibrium described by Eqs. (4.7), in which all quantities vary slowly with time. The fluid, however. 


is not taken to be in a strict hydrostatic equilibrium, which would imply the complete absence of internal motions. 
Instead, the fluid is assumed to be in an irrotational state, which implies the existence of a velocity field described by 
Eq. (4.4). The internal motions are established even when the time dependence of the tidal field is arbitrarily slow. 
By contrast, the strict hydrostatic equilibrium previously studied by Damour and Nagar [28] and Binnington and 
Poisson can only be established when the tidal field is strictly time-independent. As such, it represents a much 
less realistic configuration for the perturbed fluid. 


V. GRAVITOMAGNETIC LOVE NUMBERS OF AN IRROTATIONAL COMPAGT BODY 


In this section we determine the impact of the internal motions described by Eq. (4.4) on the body’s response to an 


applied tidal field, as measured by the gravitomagnetic Love numbers A™’’'®; we shall show that they have no impact 
on the body’s gravitoelectric response. We assume that the tidal environment varies slowly with time, and neglect all 
time derivatives in the equations that determine the metric perturbation Pap. But while the slow time dependence 
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is assumed not to have an impact on the field equations, it is crucial in the establishment of the internal motions, as 
was explained in the preceding section. 

To proceed we introduce a decomposition of into tensorial spherical harmonics. Denoting := (t, r) and 
9^ = (0, (j)), we have 


Im 

PaB = Y. + Y 

im im 

PAB = T^Yi^'^'^ABY^'" + + E 


im 


im 


in which hi^, hi^, and hi^ depend on only, and 


Fj"* := DaY^”^, 


\rim . _ 

^ AB • — 


DaDb^A^I^Y)^ab 


Y 


im 


-\rtm . ^ Jd Ti 'w^m 

Aa ■= -£a ^b^ , 

Xi^:=^(DAXi"^ + DBXiA 


(5.1a) 

(5.1b) 

(5.1c) 

(5.2a) 
(5.2b) 
(5.2c) 
(5.2d) 


are the vectorial and tensorial harmonics constructed from the standard spherical harmonics F^™(6*"^); Oab ■= 
diag(l,sin^ 9) is the metric on a unit two-sphere. Da is the covariant derivative operator compatible with this metric, 
£ab is the Levi-Civita tensor on the unit two-sphere, with components = sin0, and u pper- case Latin 

indices are raised with 11'^^, the matrix inverse to . The terms involving andF^’g inEq. (5.11 constitute 

the even-parity sector of the perturbation, and the terms involving and X^ constitute the odd-parity sector. 

The body’s response to the applied tidal field is measured by the Love numbers Kf and introduced by 

Binnington and Poisson |29j . Because these are gauge-invariant in the usual sense of perturbation theory, we may 


Aim 


calculate them in any gauge, and for this purpose it is convenient to adopt the Regge-Wheeler gauge, for which 
G^"*, and / 12 ™ are all set equal to zero. The field equations further imply that /if™ = = 0, and we find that 

Eqs. (4.4) become 


= 0 , 


VA 


= -Xe-^ptA = -Xe-'l’ Y 


'■im 


(5.3) 


im 


in Regge-Wheeler gauge. This implies that the internal motions associated with the irrotational state affect only the 
odd-parity sector of the perturbation. This, in turn, implies that the body’s gravitoelectric response to the tidal field 
is unaffected by the internal motions, but that there is an impact on the gravitomagnetic response. In other words, 
the internal motions do influence the gravitomagnetic Love numbers but they leave the gravitoelectric Love 

numbers unchanged with respect to a strict hydrostatic equilibrium. We have inserted a factor A = 1 in Eq. (5.3) 
to track the impact of the internal motions on the computation of setting A = 0 would place the fluid in a 

strict hydrostatic equilibrium. 

The computation of the gravitomagnetic Love numbers proceeds as detailed in Ref. [53] and streamlined in Ref. [33| . 
The only change concerns the additional contribution to STap, the perturbation of the fluid’s energy-momentum tensor, 
created by the internal motions described by Eq. (5.3). Focusing on the odd-parity sector, writing 


STcp = ip+ p){UaSup + UpSUa) + PPaP, 


and making use of Eqs. (3.1) and (5.3), we find that the only relevant component is 


6TtA=[{X-l)p + Xp]Yhi^X‘A^- 


(5.4) 


(5.5) 


im 


With this, the Einstein field equations linearized about the unperturbed solution of Sec. [ll| imply that the perturbation 
fields /if”*(r) are determined by 


2^2/if™ d/if™ 


Qhf^ = 0, 


(5.6) 
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FIG. 2. Rescaled gravitomagnetic Love number for £ = 3. The Love numbers are computed for selected polytropes with 
polytropic index n, and plotted as a function of 2M/R, up to (and slightly beyond when n > 1) the maximum value for the 
given equation of state. The positive values correspond to A = 0, the strict hydrostatic equilibrium. The negative values 
correspond to A = 1, the irrotational state. The circled data points at 2M/R = 0 are the post-Newtonian values obtained in 
Sec.r^ 


where 


P:=47rr^/ 

Q := f~^ [£(£ + 1) — 4m/r — (2A — l)87rr^(/r -h p)]. 


(5.7a) 

(5.7b) 


We observe that the switch from strict hydrostatic equilibrium to an irrotational state changes the sign of the 87rr^(p-|- 
p) term in Q. 


solution is matched at r = i? to the external solution provided by Binnington and Poisson 


Equation (5.6) is to be integrated outward from r = 0, near which the solution behaves as oc The internal 

(5.8) 


hi^ = 




3(£- 1)£ 


Aa-2 


£ + 1 




where 


A 4 ■= F(-£ + !,-£- 2; -2£; 2M/r), 
B 4 := F(£- 1,£ + 2; 2£ + 2; 2M/r) 


(5.9a) 

(5.9b) 


with F[a,b]C]x) denoting the hypergeometric function, and is the spherical-harmonic packaging of the tidal 


moments B 




introduced in Eq. (1.11, defined by 




£m 


(5.10) 




























FIG. 3. Rescaled gravitomagnetic Love number for = 4. The caption of Fig. [^provides additional details. 


in which := [sin 9 cos cj), sin 6 sin (f), cos 9]. 

A practical formulation of this procedure was provided by Landry and Poisson [33]. It involves introducing the 
logarithmic derivative K( := d\nhf^/dlnr, which satisfies the nonlinear differential equation 

r-^ + Ki^Ki — 1) — Pk^ — (5 = 0. (5.11) 

This equation is integrated from r = 0, at which = .^ + 1, to r = i?, at which ki = k^. The matching condition at 
r = R produces 


n \2t 


with 


r mag _ 

hj/) — 


e R RA'^-{K\-e- 1)^4 
2 {i + 1) 2M - [ki + e)Bi ’ 


(5.12) 


(5.13) 


in which a prime indicates differentiation with respect to r. The overall scaling of with (i?/2M)^^ differs 

from the convention adopted in Damour and Nagar [28] and Binnington and Poisson [29] . who introduced instead a 
scalefree Love number via the relation AT™'^® = . The missing factor of R/2M was incorporated 

in Eq. (5.13), and it ensures that fc“^® tends to a nonzero value in the limit 2M/R —>■ 0; this is unlike which 

goes to zero in the limit. 

In Figs. [l]|^ we present sample calculations of the rescaled gravitomagnetic Love numbers for a polytropic 
equation of state p = Arp^+^Z”, e = np, where AT and n are constants. The details of integrating Eq. (5.11) for this 
specific case are presented in Appendix [^ We see from the hgures that the tidal response of an irrotational body 
is dramatically different from the response of a body in strict hydrostatic equilibrium: the gravitomagnetic Love 
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FIG. 4. Rescaled gravitomagnetic Love number for £ = 5. The caption of Fig. [^provides additional details. 


numbers of an irrotational body are negative, while they are positive for a body in strict hydrostatic equilibrium. The 
change of sign can be seen to originate from the 87rr^(/r + p) term in Q. Both sets of Love numbers, however, share 
the properties that they decrease (in absolute value) with increasing 2M/R and with increasing polytropic index n; 
because increasing n decreases p for a given p, and therefore produces a body that is more centrally dense, both 
properties are associated with the fact that a more compact body develops smaller multipole moments. 


VI. POST-NEWTONIAN CALCULATION OF THE GRAVITOMAGNETIC LOVE NUMBERS 


In this section we exploit the methods of post-Newtonian theory to calculate the gravitomagnetic response of a 
weakly self-gravitating body to an applied tidal field. Our calculation returns an expression for that reproduces 

the leading-order behavior of the relativistic function when the body’s compactness 2GM/{c^R) is small. In this 
section we restore powers of G and c, which were set equal to unity in the preceding sections. 

The post-Newtonian metric is expressed in harmonic coordinates as 

goo =—^ + -^U + 0{c ^), goj = --^Uj + 0{c ^), gjk = + 0{c ^), ( 6 - 1 ) 


and it involves a vector potential Uj in addition to the familiar Newtonian potential U , which satisfies the Poisson 
equation = —47rGp. The post-Newtonian terms of order c~'^ in goo are not required in our developments, and 
the second post-Newtonian terms of order c~^ in goj are absorbed into Uj. 

As was described in Sec.|Tj the tidal environment is taken to be a pure gravitomagnetic Apole field characterized by 
the symmetric-tracefree (STF) tensor Bk.^k 2 ---ki{t), which is assumed to vary slowly with time, so that we may neglect 
time derivatives in the field equations. The tidal field is described by the vector potential 


^tidal _ _ 


6(€- I) 


r e 


JPq^^^‘^k2k3. 




( 6 . 2 ) 
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where ejpq is the completely antisymmetric permutation symbol, := SjkX^x^, and 51-1 := x^/r. The STF nature of 
the tidal moment Bkik 2 ---ke{t) ensures that the vector potential satisfies Laplace’s equation = 0 as well as the 

harmonic gauge condition = 0. It also implies that the vector potential can be expressed in the alternative 

form 


^tidal _ _ 


6 (^- 1 ) 


feafea- 


■kt 


^{pk2k3---kf) 


(6.3) 


in which 57 * 1 ^ 7^2 ... ^kt angular brackets indicate the operation of trace-removal. It is ea sy t o 

show that a transformation to spherical polar coordinates (r, 9"^) brings the metric to the form displayed in Eq. (5.8|, 
in which we set A 4 = 1 + 0(c“^) to respect the post-Newtonian approximation, and B 4 = 0 because the body’s 
response has not yet been incorporated within the vector potential. 

The body’s response to the applied tidal field is governed by the post-Newtonian Euler equation and the field 
equation satisfied by the vector potential. The relevant aspects can be collected from the textbook by Poisson and 
Will [55]. From their Eq. (7.15b), in which we neglect the time derivatives, and their Eq. (7.23b) we obtain 


U, = 


tidal 


G 




(6.4) 


in which the domain of integration ^ is truncated to r' < TZ, where 7^ is a cutoff radius to be specified below. 
The body’s response is captured by the Poisson integral, and the integrand involves the effective current density 
s-l := which incorporates both a matter and field contribution. An expression for this is displayed in Exercise 

8.4 of Poisson and Will; we have 


= pv^ [1 + 0(c-2)] + - d^W)dkU + 0(c-4), (6.5) 

in which we neglect a term involving dtU in the field contribution to the mass current, as well as the 0 {c~^) terms 
multiplying pv^ because, as we shall see presently, the velocity field is itself of order c“^. Our expression for Uj should 
be multiplied by the factor (1 — 2Ujc?) arising from Eq. (7.23b) of Poisson and Will. This factor can be neglected 
when it multiplies the Poisson integral, but it produces a relativistic correction to We shall not be interested 

in such corrections, and choose to drop the factor all together. 

The velocity field is determined by the irrotational condition that we wish to place on the fluid. Writing = 7 ( 0 , v^') 
for the four-velocity, in which 7 is a normalization factor, a short computation reveals that the circulation of a post- 
Newtonian fluid around any spatial circuit c is given by 


C{c) = 'jh 

An irrotational state therefore requires 


1 -I- -^U \ vj — -^Uj + 0{c ) 


, = ^C/,+0(c-4). 


dx^. 


( 6 . 6 ) 


(6.7) 


This expression agrees with results previously obtained by Shapiro |30j and Favata m- 


The calculation of the body’s response proceeds by inserting Eq. (6.71 within Eq. (6.51, and replacing Uj by 


jjtidai _|_ 2 ^ jg expression. We find that the velocity contribution to the mass current is given by 


»<? r){pk2k3-kt) 

pi tjpqU k2k3--kt‘‘ > 


and that the field contribution is 


2A 

' ~3(£- l)c2 

.._ ^ ^ nipUK3--- 

67r(£-l)c2”' 


kf) 


( 6 . 8 ) 


(6.9) 


To arrive at the last expression we recalled the fact that for a spherical body, dkU = —{Gm/r^)ilk, in which m(r) is 
the mass inside a sphere of radius r. We also inserted a factor of A = 1 in Eq. ( 6 . 8 ) to again track the influence of 
the internal motions on the body’s gravitomagnetic response. 

The next step is to evaluate the Poisson integral in Eq. (6.4). The calculation requires the addition theorem for 
spherical harmonics. 


1 


\x — x'\ 


4_ 

£m > 


( 6 . 10 ) 
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in which r< = min(r, r') and r> = niax(r, r'), and the identity displayed in Eq. (1.171) of Poisson and Will, 

^ Y,m{9, f yi^{9\ dY)! = (6.11) 

m 

Evaluation of the Poisson integral for sj is straightforward, and when r > R we obtain 




_ I _OPR? Qk2Qk3 

0 ('£_X) “ 




16A ^ GM 


2 £ + r 


( 6 . 12 ) 


in which 


Mi? 2 ^ 


pr 


2t+2 


dr 


(6.13) 


is a dimensionless moment of the mass density. In this case the domain of integration is naturally limited to the 
volume occupied by the body, and the cutoff radius TZ is irrelevant. It does, however, play a role in the Poisson 
integral for s|, which features the radial integral 


rU 


J := 


m{r')r 






dr'. 


(6.14) 


To evaluate this for r > R we break up the integration domain into a first segment 0 < r' < i? in which 

r> = r, and m(r') is unspecified, a second segment i? < r' < r in which r< = r', r> = r, and m{r') = M, and a third 

segment r < r' < TZ in which r< = r, r^ = r', and m{r') = M. We obtain 


J = 


1 


.1+1 


[ m{r)r‘^‘^-^ dr -—MR^^ 

Jo 


2^+1 M ^ 

H- Mr" " - r" 

2i TZ 


(6.15) 


To simplify this we integrate the first term by parts, making use of the Newtonian field equation dm/dr = Airr'^p. 
We also discard the term proportional t o M r"~^, because it merely gives rise to an uninteresting correction of order 
GM/{c^r) to the tidal potential of Eq. (6.2). The term proportional to Mr"/TZ can be seen to alter the amplitude 
of the tidal potential by a correction of order GM/{c^TZ), and we eliminate this meaningless shift by setting TZ = oo. 
With all this we find that 


J = - 


2 ^ MR 


21 


rt+l 


and completing the evaluation of the Poisson integral, we arrive at 


u^ = - 


1 




0 ('£_ X) “ 




8 (^ +1) GM R^" 1 


(6.16) 


(6.17) 


Adding the contributions, we find that the body’s gravitomagnetic response to the applied tidal field is described by 


rrresp __1_ I nPnl Qk^nki 

- 0('^ _ X) " fc2fc3-fcf“ “ 




8[(2A - 1)£ - 1] GM R 


21 


42 £+l) 


(6.18) 


The complete vector potential is , and comparison with Eq. (1.1) reveals the post-Newtonian expression 


for the gravitomagnetic Love numbers. Recalling Eq. (6.13), we find that 

21 


j^mag _ I R 
" \2GM/c^ 


jLinag 


with 


f.mag_ 

t^D ■ — 


2[(2A-1)£-1] ^ r 

{£+l){2e + l) MR^" Jo 


pr 


, 2£+2 


dr. 


(6.19) 


( 6 . 20 ) 


We recall that A keeps track of the internal motions associated with the fluid’s irrotational state. Setting A = 1 places 
the body in the irrotational state, and we observe that negative gravitomagnetic Love numbers must be assigned to 
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such a body, irrespective of the multipole order £ > 2 and the equation of state. By contrast, setting A = 0 places 
the body in a strict hydrostatic equilibrium, and such a body necessarily comes with positive gravitomagnetic Love 
numbers. These properties were featured in the fully relativistic results displayed in Figs. [T]|^ and indeed, we observe 
that the post-Newtonian expression of Eq. (6.20) accurately reproduces the relativistic Love numbers in the limit 
2GMI{c^R) ^ 0. 

Our post-Newtonian calculation of the gravitomagnetic Love numbers completes previous attempts carried out by 
Favata m and Damour and Nagar [^. In his work (see Sec. Ill B of his paper), Favata introduces a definition for 
the Love numbers that accounts only for the velocity term s) in the effective mass current; it omits the field term s^. 


and Favata therefore produces only the A term in Eq. (6.20). On the other hand, the post-Newtonian calculation of 


Damour and Nagar (see Sec. VIII of their paper) places the fluid in a strict hydro static equilibrium instead of the 
irrotational state, and therefore accounts only for the A-independent term in Eq. ( 6.20| ). (It should be noted that 
the Damour-Nagar definition for the gravitomagnetic Love numbers includes a multiplicative minus sign compared 
to ours; their Love numbers are negative when ours are positive.) Our own calculation brings the two partial stories 
together, and generalizes the previous calculations (which were limited to i = 2) to arbitrary multipole order i. 

Equation (6.20) can be evaluated in closed form in a few simple cases. For a constant density body we find that 


rmag _ _ 


3[(2A- 1)£- 1] 


2{e + i){2e + i){2e + 3 )' 

For an n = I Newtonian polytrope, for which p = M/{AE?r) sin( 7 rr/i?), we have that 

~ 207r" + 120 ) _ (_0.0068498, 0.020549), 
zOtt^ 


7 mag _ _ 
^3 “ 


7 mag _ 
— 


^mag _ _ 


~ (-0.0044829,0.0089658), 


(3A - 2)(7r6 - 427r4 -h 8407r2 - 5040) 

287r® 

(8A - 5)(7r8 - 727r6 -h 30247r4 - 604807r2 -h 362880) 

OOtt® 

(5A - 3)(7rio - IIOtt® -h 79207r6 - 3326407r4 -h 66528007r2 - 39916800) 


~ (-0.0028102,0.0046836), 


( 6 . 21 ) 

( 6 . 22 a) 

( 6 . 22 b) 

( 6 . 22 c) 


6671^® 

The numerical values correspond to A = 1 and A = 0, respectively. 


-(-0.0018302,0.0027454). (6.22d) 
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Appendix A: Structure and perturbation equations for polytropes 


To perform the computations described in Sec. for the polytropic equation of state p = e = np we 

recast the background field equations of Sec. [IT] and Eq. (5.11) in convenient, dimensionless forms. To achieve this 
we introduce the central density pc ■= p{r = the central pressure Pc ■= p{r = 0), the length scale vq defined by 
Tq := (n -I- l)pc/( 47 rp^), and the mass scale mg := 47rpc?'o- A useful dimensionless parameter is b := pdPc = Kp]J^^ 
which can act as a substitute for the central density as a label of polytropic models. A frequently encountered 
combination of scaling quantities is mo/rg = {n + 1 ) 6 . 

We next introduce the dimensionless, Lane-Emden-type variables 0, and iz, such that r = rg^, p = 
p = Pc0"'~^^, p, = Pc0"(l + nb9), and m = mg^^zz. The equations that determine the internal structure of the 
unperturbed polytrope are then 


^ = r (1-f n66<) - 3zz 


(Al) 


and 


?^ = -er^ [1 + (n + 1)60] (zz + 60"+i), 


(A2) 
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with / = 1 — 2(n+ An equation can also be displayed for the gravitational potential ip, but this is not needed 

to calculate the Love numbers. The integrations begin at ^ = 0 with 0(0) = 1 and 1 /(0) = |(1 + nh). They proceed 
until ^ at which 0 = 0 and v = Vs. The body’s compactness can then be calculated as 2M/R = 2{n + z/g- 
In the limit 6 —>■ 0 the equations reduce to the standard Lane-Emden form. 

The dimensionless version of Eq. (5.11 1 is 


^— + - 1) - Pne - Q = 0, 


(A3) 


with 


P={n + 1)6^20"/-! [1 + (n + 1)60], 

Q = i |f(£ + 1 ) - 2(n + 1)6^" [ 21 / + (2A - l)0"(l + {n + 1)60) 


(A4a) 

(A4b) 


The integration begins at ^ = 0 with Kg = £ + 1 and proceeds until ^ at which 

In practice it is helpful to use x := In^ as the independent variable, and to start the integration at a large, negative 
value of X. Starting values for 0, v and can be obtained from the Taylor expansions 0 = 1 + 02 ^^ + 04 ^^ + • • •, 
v = |(1 + nh) + + • • •, and = £ + \ + + • • •) in which the various coefficients can be 

determined from the differential equations. 
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